GENERALIZED STIRLING TRANSFORM 
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Abstract. In this paper, algorithms are developed for computing the 
Stirling transform and the inverse Stirling transform; specifically, we 
investigate a class of sequences satisfying a two-term recurrence. We 
derive a general identity which generalizes the usual Stirling transform 
and investigate the corresponding generating functions also. In addi- 
tion, some interesting consequences of these results related to classical 
sequences like Fibonacci, Bernoulli and the numbers of derangements 
have been derived. 



1. Introduction 

The Stirling numbers arise frequently in mathematics, especially in enu- 
merative problems. This is the reason of their important role in combina- 
torial analysis, number theory, probability, graph theory, calculus of finite 
differences and interpolation. The notations for these numbers have never 
been standardized, this paper follows the notation of Riordan for the signed 
Stirling numbers of the first kind s (n, k) and Knuth's notation for the Stir- 
ling numbers of the second kind {^}. 

The Stirling transform of a sequence (o n ) is the the sequence (6 n ) given 

by 

(1.1) b n = f:{i}a k , 

and the inverse transform is 

n 

(1.2) a n = (n, k) b k . 

k=0 

The identity (jl.ip has a combinatorial interpretation given in [2]. If a n 
is the number of objects in some class with points labeled 1,2, ... ,n (with 
all labels distinct) then b n is the number of objects with points labeled 
1, 2, . . . , n (with repetitions allowed). 

In this paper, algorithms are developed for computing the Stirling trans- 
form and the inverse Stirling transform; specifically, we investigate a class 
of sequences satisfying a two-term recurrence. We derive a general identity 
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which generalizes the usual Stirling transform and investigate the corre- 
sponding generating functions also. 

Given a sequence a m := ao. m (m > 0). We construct an infinite matrix 
S := (a n>m ) as follows: 

The first row ao, m of the matrix is the initial sequence; the first column 
b n := a n> o (n > 0) is called the final sequence and, each entry a„ )m is given 
recursively by 

(1-3) Cl n -)_i )m = CLn,rn+l ~\~ ^ a n,m- 

Conversely, if we start with the final sequence, the matrix S can be re- 
covered by the recursive relations 

(1-4) a n ,m+l = CLn+l,m ~ m&n,m- 

2. Definitions and notation 

In this section, we introduce some definitions and notations which are 
useful in the rest of the paper. N being the set of positive integers and 
N = NU {0} . 

The falling and rising factorials are defined, respectively by 
(x) n = x (x - 1) • • • (x — n + 1) , (x) = 1 

and 

(x) n = x {x + 1) • • • (x + n — 1) , (x} = 1. 

The (signed) Stirling numbers s(n,k) of the first kind, which are usually 
defined by 

n 

(2.1) (x) n = Y,s(n,k)x\ 

k=0 

or by the following generating function 

(2.2) i f (In(l + x))* = ^*(n,fc)^ r . 

n>k 

It follows from (j!TTj) or ([272]) that 

(2.3) s (n + 1, k) = s (n, k — 1) — n s (n, k) 
and that 

s(n, 0) = 5 n fi (n € N) , s (n, fe) = (k > n or k < 0), 

where (5 n)m denotes the Kronecker symbol. 

The Stirling numbers {^} of the second kind count the number of possible 
partitions of a set of n objects into k disjoint blocks. These numbers can be 
defined explicitly by 



■=e{;}w. 

k=o ( J 
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For any positive r S N the quantity {)!} r denotes the number of partitions 
of a set of n objects into exactly k nonempty, disjoint subsets, such that the 
first r elements are in distinct subsets. These numbers obey the recurrence 
relation 



0, n < r, 



r 



(2.1) {"J =6 ktr , n = r, 

n] fn — 1") [n — 1] 



The exponential generating function is given by 



(2.5) El? + r i - = \-/ x {e x -\f 

n>k ' 



The properties 



r n-r 



and 



r-l 



are given in [3], which one can consult for more details on r-Stirling numbers. 



3. Combinatorial identities 

Theorem 1. Given an initial sequence («o,m) m >o , define the matrix S by 
111. Sty . Then, the entries of the infinite matrix S are given by 



(3.1) 
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Proof. We proof by induction on n, the result clearly holds for n = 0. By 
induction hypothesis 



E (n + m + 1) (n + m) 

s , , ao, m +k+i + m< > 

i—n I K + 771 + 1 J m+1 I m J I 



El n + to 



Ein + m + l) (n + m 

i , . . -, f a , m +fc+i + to<^ 
I A; + m + 1 J +1 [ to 



fc=o 

n-2 
k=Q ^ 



ao,: 



n + m 
k + m + 1 



n + to + 1| (n + m + 1') 

n + m+1) m+1 i « + m J m+ i 

(n + m) f n + m ] 

+ m< ) ao^m + m) < > a 0m+h+ x. 

I m )m ^\ k + m + 1 )m 



From (|2.6p and after some rearrangements, we get 



n+l 



n + m+1 



0"n,m+l + fno^n^m — j. , f a 0,m+A:- 

fc=0 11 



□ 



Theorem 2. Given a final sequence (a n ,o) n >o , define the matrix S by jl.4\ )- 
Then, the entries of the infinite matrix S are given by 



in 

(3.2) a„, m = (to, k) a n+k:0 . 

k=0 
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Proof. We proof by induction on m, the result clearly holds for n = 0. By 
induction hypothesis and (|2.3p . we have 

m+1 m 

a n +i,m - ma n , m = ^ s (m, k - 1) a n+ k,o - rn^s (m, k) a n+ k,o 



k=l k=0 

in 

■■ s (m, m) a„ +m+ i,o + ^/ (m, k - 1) a n+kfi 
k=l 

in 

ms (m, 0) a nfi - m^s (m, fc) a n +fc,o 

k=l 

m 

■■ s (m, m) a n+m+ i j0 + ^ (s (m, fc - 1) - ms (m, fc)) a n+fc)0 
fc=i 

ms (m, 0) a n) o 

□ 



Corollary 1. 

(3.3) 



Z/( m > fc ) 6 «+* = X)u 4. Z\ am+k - 



The last identity can be viewed as the generalized Stirling transform which 
reduced, for m = 0, to the Stirling transform (jl.ip of the sequence a n , and 
for n = reduces to the inverse Stirling transform (jl.2p of the sequence b m . 
We may now formulate the following algorithms 

Algorithm 1. Stirling transform 

Input: a n 

Output: b n 

Set X m — a, n — Ta m 

for n = 0, 1, . . . do 

for m = n, n — 1, . . . , do 
X m -\ := (m — 1) X m _i + X m 
end do 

end do 

Algorithm 2. inverse Stirling transform 
Input: b m 
Output: a m 

Set Y n — b n ^ m — n 
for m = 0, 1, . . . do 
Y — b 

for n = m, m — 1, . . . , do 
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Y n -i := Y n -(m- n) F n _i 
end do 
a m := Y 
end do 

Example 1. Setting ao, m = 1 in ( fff.ffj) . we get the well known identity [8] 



E 

fc=0 



n + m 
k + m 



y~]s(m, k)B n+k , 



k=0 



where B„ is the nth Bell number. 



Example 2. Let (F n ) nl - No be the Fibonacci sequence given by Binet's for- 
mula 

where a = 1+ 2 and (3 = 1 ~ 2 V ^ , If the initial sequence ao,m = 
then we get the following matrix 



«-«> m - HO J 





i 





1 


-4 


19 


-108 


1 


i 


1 


-1 


3 


-13 


71 


1 


2 


1 





-1 


6 


-37 


2 


3 


2 


-1 


2 


-7 


34 


3 


5 


3 


-1 


1 


-1 


-3 


5 


8 


5 


-2 


3 


-8 


31 


8 


13 


8 


-3 


4 


-9 


28 


13 


21 


13 


-5 


7 


-17 


59 


V: 















/ 

From this matrix we observe that a n o = a n ^ = — fln+2,3 = F n , and 
= L n , where (-^n) ng N the Lucas sequence given by Binet's formula 

L n = a n + f5 n . 

It is well known that the F n and L n are connected by the formula 

L n = F n ^x + F n+1 , (neN). 
By h3. 3\) . one can deduce that 



and by Theorem^ we get 

2 

Fn = J] S ( 2 > k ) F n+k = -F n +1 + F n+2 



■a, 



m+k 



( P)m+k) ' 



k=0 



-/^ (3, k) F n+2 +k — —2^+3 + 3F n+ 4 — Fj 



n+5) 



k=0 
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and for n G No 

4 

L n = ^/ (4, fc) F n+3+k = -6F n+4 + HF n+5 - 6F n+e + F n+7 . 
fc=o 

By TheoremUl 

n 



{: : 2 2 } 2 



and 



fc=o 

n+3 



4. Generating function 

Theorem 3. Suppose that the initial sequence ao, m +r has the following expo- 
nential generating function A r (z) = ^^CLo,k+r^r- Then the sequence {a n , r } n 

of the rth columns of the matrix S has an exponential generating function 
Br (z) = y^a n r — - given by 

n>0 

(4.1) B r (z) = e rz A r {e z - 1) 



Proof. We have 



B r (z)=J2*0,r +k j:{ n k l r r } ^ 

fc>0 n>0 k ' -> r 

= Yj a °>r+ky } e rz {e z - 1) 



fc>0 



(e z - 1) 



+fc " fc! 

fc>0 



e r2 A r (e z - 1) 



□ 



Theorem 4. Suppose that the final sequence a n+r fl has the following expo- 
nential generating function B r (z) = ^^ak+rfijr ■ Then the sequence {a rtm } m 

k>0 



8 



MOURAD RAHMANI 



of the rth rows of the matrix S has an exponential generating function 



m>0 



given by 



(4.2) 

Proof. We have 



A r (z) = B r (lll(l + Z)) . 



k>0 



m>0 



z m 

A {z) = ^a r+fci0 ^2 s ( m ' k ) ~ 



(ln(l + z)f 



y~^Qr+fc,0 
fc>0 

B r (ln(l + z)). 



A:! 



□ 



Example 3. A derangement on a set {1,2,... ,m} is a permutation ir = 
*x*2 ■ ■ ■ *'m such that i k ^ k for k = 1, 2, . . . m. The number of derangements 
on {1,2, . . . ,m} is denoted by D m and given by D m = [— 1 , where [x] the 
nearest integer function. Now, if the initial sequence ao, m = (— l) m D m , then 
we get the following matrix 

•A 



(1 





1 


-2 


9 


-44 





1 





3 


-8 


45 


1 


1 


3 


1 


13 


-39 


1 


4 


7 


16 


13 


76 


4 


11 


30 


61 


128 


159 


11 


41 


121 


311 


671 


1381 


V 













/ 



The generating function of the sequence ao, m is Aq(z) = It follows from 
\4--l\ l that Bq(z) = exp(e z — z — 1), and we notice that a n] o is the number v n 
of partitions of {1,2, ... , n} without singletons (see for instance [10 J. By 
\3. 3\) . one can then deduce that 

m n , v 

^2s(m,k)v n+k = Y J (- l ) m+k W Z\ D m+k 
k=0 k=0 ^ > m 



k=0 

If m = 0, we have 



k=0 k ^ 



Example 4. T/ie exponential generating function of the Bernoulli polyno- 
mials B n [x) is 



So (z) := = J> (x) 



n>0 
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By Theorem 4, we have 

A (z) 

It is not difficult to show that 
[z m ] A (z) 



(l + zfln(l + z) 



E(-D 



i=0 



m — i + 1 ' 



where [z n ] f (z) denote the operation of extracting the coefficient of z n in 
the formal power series f (z) = ~^2f n z n . Now, let us consider S defined by 
(TOj with the final sequence a n $ = B n (x) , by h3. 3\) . we have 



(m, k) B n+k (x 

k=0 k=0 



n r m+k 

Y\l ) Eh 



m+k—i 



m + k — i + 1 



Example 5. Catalan and Motzkin numbers naturally appear in a large num- 
ber of combinatorial objects. It is well known that the Catalan number 

\n/2\ 

Cn = ( 2 ™) an d Motzkin number M n = ^rpy (^.) ( 2 k ) are connected 

k=0 

by® 



a 



n+l 



1,—n \ / i—n \ / 



fc=0 v ' fc=0 

Using the generalized Stirling transform, we can show that the Catalan num- 
bers are related with Motzkin numbers in terms of Stirling numbers by 



n+l 



(4.3) 

and 
(4.4) 



(n, k) M k = ( n + 1> fc ) 



fc=0 



fc=0 



C„ = <5„, + ™ | s (k - l,i) Mi^M n = l\ s ( k + !•*) C - 

fe=li=0 ^ ' fc=0i=0 ^ ' 

Setting the final sequence a n> o = C n , we get the following matrix 

\ 



n k+l 



( 1 


l 


1 


1 





1 


-5 


29 


1 


2 


3 


3 


1 





-1 


7 


2 


5 


9 


10 


4 


-1 


1 


-1 


5 


14 


28 


34 


15 


-4 


5 


-11 


14 


42 


90 


117 


56 


-15 


19 


-42 


42 


132 


297 


407 


209 


-56 


72 


-160 


132 


429 


1001 


1430 


780 


-208 


272 


-614 


V ! 
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Since 



B (z) 



n>0 



111 



1/2 



Az , 



where \F\ ( ; z 



^^^yz™. It follows from gjf i/iat 



n>0 



A(^) = ^i?n^=iFi( X f ;41n(l + z 

n>0 n ' ^ 



and 



(4.5) 



J2 S ( m > k ) C n+k=J2\ k + Z\ Rm+k 
k=0 k=0 ^ J m 



Now, if the initial sequence ao, m = Rm+i, we get the following matrix 





/I 


1 


1 





1 


-5 


29 


-196 




1 


2 


2 


1 


-1 


4 


-22 


146 




2 


4 


5 


2 





-2 


14 


-100 


T = 


4 


9 


12 


6 


-2 


4 


-16 


93 


9 


21 


30 


16 


-4 


4 


-3 


-26 




21 


51 


76 


44 


-12 


17 


-44 


172 




51 


127 


196 


120 


-31 


41 


-92 


282 




\: 
















this matrix 


we 06 


serve 




On,0 = 




We 


prove 



using generating functions. We have 



git 



n>0 
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From Ity , we get 

Bq(z) 




n! 



Cfc+i — r 



n! 



n>0 



follows 
(4.6) 



E s ( m > fc ) M «+ fc = Eu + Z r • 

k=o k =o l +m)m 



Combining results [4-®\) an d gives \4-3\ ) and h4-4\ )- 

5. Hankel transform 

The Hankel transform of a sequence a n is the sequence of Hankel deter- 
minants det (ai + j) Q<i - <n . A number of methods for computing the Hankel 
determinants have been widely investigated [U EJ [9] . It is well known that 
the Hankel transform of sequences a n and (5 n are equal under the binomial 
transform 17 



fc=0 v 7 



A natural question arises: "What about the Hankel transform of the 
sequences a n and b n under the Stirling transform?" In this section we show 
that there is a connection between the generalized Stirling transform and 
the Hankel determinants. 



Theorem 5. For n > 0, we have 



det (a 



Proof. We can write 
det(6 i+i ) <. i < n 



i,3JQ<i,j<n 



«o,o a i,o 

Ol.O a 2,0 



det (bi+j) 



O2,0 
^3,0 



0<i,j<n ' 



03,0 
04,0 



«n,0 On+1,0 On+2,0 «n+3,0 

after applying (11. 4p . the determinant is unchanged 



On,0 
«n+l,0 

O2n,0 
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^0,0 
Ol,0 



Using (|3.2|) . we get 

det(6i+j)o<<j-<„ 



ai,o 

O2,0 



02,0 — «1,0 
«3,0 — «2,0 



a 3,o — 02,0 — 2 (02,0 — ai,o) 

«4,0 — «3,0 — 2 (03,0 — 02,0) 



«n,0 fln+1,0 O n+ 2,0 ~ «n+l,0 On+3,0 _ «ji+1,0 _ 2 (a n+ 2,0 ~ «n+l,o) 



^0,0 Oo,l Oo,2 «0,3 • • • O0,n 

ai,o «i,2 «i,3 • • • a i,n 



dn,0 a n,l «n,2 On,3 a n 

from which the relation follows. 

The answer to the previous question is given in the following 
Corollary 2. For n £ No, we have 



□ 



det 



i +3J0<i,j<n 



det E 

\fe=0 
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